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SUMMARY 

Rosen’s treatment of impact ionization of neutral atoms is applied to the 
collision of ground- state hydrogen atoms near the Ionization threshold. The 
singlet interaction potential is approximated by a Morse potential and the trip- 
let interaction potential is approximated by an exponential repulsive potential. 
The energy- dependent Ionization cross section is presented for collision energy 
between threshold and twice threshold. An implication of the classical path 
assumption is briefly discussed. 


INTRODUCTION 

Several years ago an approximate method was introduced by Rosen (ref. l) 
for the calculation of atom-atom impact ionization cross sections for collision 
energies near threshold. At that time it was applied to helium- helium colli- 
sions (ref. l) and argon-argon collisions (ref. 2). The method is used herein 
to compute the energy- dependent cross section for the interaction 

H( Is ) + H( Is ) -> H(ls) + H + + e" (l) 


near threshold. 

This calculation has a similar but simpler perturbation matrix element 
than those considered by Rosen; hence, much of the calculation required to ob- 
tain a cross section for process (l) has already been reported (ref. l). 

The only previous calculations (ref. 3) for process (l) made use of the 
Born approximation and therefore only apply at energies considerably above 
threshold. 


ROSEN’S APPROXIMATION 

The method assumes that the two colliding atoms follow classical paths 
dictated by their interaction potential, but at a time corresponding to their 



distance of closest approach a perturbation is suddenly turned on. The pertur- 
bation matrix element Vf ^ between the initial state i and final state f 
of the two-atom system decays very slowly because the incident atom has lost 
most of its kinetic energy and is moving away very slowly. Obviously such con- 
ditions will be met only if the kinetic energy of the incident atom is no more 
than about twice the threshold value. The result of an ionization cross-section 
calculation under these conditions is obtained from reference 2 as 


a 



b db 


( 2 ) 


where I is the ionization energy, e = l/2 (pv^) is the relative energy of col- 
lision, m is the electron mass, p is the ejected electron 's momentum, and 
Vfi(R) is the perturbation matrix element. (All symbols are defined in the 
appendix. ) The relation between the distance of closest approach R and the 
impact parameter b is given by the classical expression 
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(3) 


where U(R) is the interaction potential between the two atoms. 


Perturbation Matrix Element 


The perturbation matrix element for the interaction 

H(ls) + H( Is ) H( Is ) + H + + e" 


is written 


V 


fi 



dT 


(i) 

( 4 ) 


where V is the perturbation potential, the product of two Is wave 

functions and is the product of a Is wave function and a plane wave 

representing the freed electron. 



The coordinate system employed in the cal- 
culation is shown in figure 1. Atom A is con- 
sidered to be incident on atom B which is at 
rest. Electron 1 is initially bound to atom A ; 
electron 2 is initially bound to atom B. The 
perturbation potential is then given by 
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V 


(5) 


e 2 + e 2 e 2 e 2 

R r 12 ^ B 1 


The initial and final wave functions are 

= cp(A 1 )cp(B 2 ) 

(ip-B 2 )/ft 

tf = cp(A 1 )e 

respectively, where cp is the normalized Is wave function given by 


( 6 ) 


„(r) = (£) 


3/2 -|_ -(z/a )] 


y/n 


(7) 


where a Q is the Bohr radius and z is the effective charge. Combining equa- 
tions (4) to (6) gives 


V^. = — 
fi B 


J cp(A 1 )cp(A 1 )cp(B 2 )e^ 1P B 2^/ fi ax n dx 


T l> t 2 


1 ua 2 


+ e 


/ 


q3(A ] _)cp(A 1 )cp(B 2 ) — e 


1 (-ip*B 2 )/ft 


12 


dx-^ dx ^ 


xp t 2 


- e 


/* l ?(A 1 )cp(A 1 )cp(B 2 ) -L e^ 1P B 2^/ fi dT _^ dT ^ 

%/ 2 


T l» t 2 


^ cp(A 1 )cp(A 1 )cp(B 2 ) ^ e (-iP*B 2 )/ ft dTi dTj . 


( 8 ) 


T l» T 2 


Taking advantage of normalization and the following integrals worked out 
by Rosen (ref. 1) 


= 4^ / cp(B 2 )e ( " ip ’ B 2 ) / tt dx. 


J 2 = e2 


/ 


cp(A 1 )cp(A 1 )cp(B 2 ) — — e 


T l’ T 2 


1 (-iP'B )/ft 

— e ^ 

r 12 


dx-^ dxg 
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/ 


cp(B 2 ) ( -ip'Bg )/fi 

H e 


dT^ 


I 4“ 


/ 


cp(A 1 )cp(A 1 ) 


B_ dT l 


give the perturhation matrix element as 

V fi = I J 1 + J 2 - X 5 - | Vt 

Rosen (ref, 1) made the approximation, Ig " I5 ” This implies that the 
atomic separation R is much greater than the electron orbital radius. Using 
this approximation and the evaluation of 1^ and 1^ in reference 1 gives 


5/2 


V fi (R) = 


8 e 2 WJLV' h + JL R \ 
VW V a o / 

fe) 1 * (# 


2(z/a 0 )R 


(9) 


Substituting equation (9) into equation (2) and defining X 2 = 

X 2 = j - 1 
m I 


a 


-&> 




2ml 

~xzl* 


{or = unity for Is hydrogen) 


r o = C-/V o ( ° 0 = 256 a o for Is hydrogen) 

° (2ml) 5 / 2 i 2 0 ° 


f(V a ) = 


/ 


X 2 dx 


(X 2 + 1) 2 (X 2 + a 2 ) 


( 10 ) 


and 


4 



g 


b db 


( 11 ) 


/ 


L- p\V 4 ( z / a o) R 


( 1 + t B ) e 


R 


yields 


a - o 0 HJ^l)e 


( 12 ) 


An integral table (ref. 4) was used to obtain the following evaluation of equa- 
tion (10): 


f(X m ,!) =- 




Xm 


7 ^n 


10 ( X m + ^ 80 ( X * + l) 480(X 2 + l)' 


7X. 


m 


7 \ 


-1. 


384( X 2 + l) 2 256(4 + lj + 256 ^ ^ 


Interaction Potential 


In order to evaluate g the variable of integration is changed from b 
to R (ref. 1). Equation (3) can he utilized to write b db as 


b db 



1 R dU(R) 

2 € dR 


R dR 


(13) 


The classical collision paths of two hydrogen atoms each in a Is elec- 
tronic state can follow one of two well-known potential energy curves corre- 
sponding to the and (parallel and antiparallel spins) states of the 

quasi-molecule (Slater, ref. 5). These potential curves are shown in figure 2 
as reported by Dalgarno and Lynn (ref. 6). 


Also shown in figure 2 are the convenient approximate forms of the poten- 
tials used for computation of g. The 3 Z u state is represented by an exponen- 
tial repulsive potential 

U 3 (R) = cce-P R (14) 

where a = 3. 6 Rydbergs and p = 1. 48 a”-*-. The -*-Eg state 
approximated by the Morse Potential (Slater, ref. 5). 

UjM = D[e- 2B < K ->-> - 2e- B < E -ri] 


is well- 

(15) 
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1 . 0 , 


V Eq. (14) 



A\ I I I I I I I I 

0 12 3 4 

Interatomic distance, R, a 0 


Figure 2. - The interaction energy of H(ls)+H(ls): solid lines (ref. 6). Dashed lines, approx- 
imate functions. 


where D =0.35 Rydberg, B = 1.08 a“\ and y = 1. 5 a c . When the poten- 

tial (eq. (14)) is used in equations (11) and (13) to compute g the result is 
designated g^; when the potential (eq. (15)) is used the result is desig- 

nated g-^. 


Evaluation of gj and g^ 

Combination of equations (11) and (13) with equations (14) and (15) gives 


/»00 

A)3 


:i + R) 2 e“ 4R 


fl - - e _(3R + i ^ Re -f3R ]dR (16a) 

L £ <~i €. 


and 


g l 


/> CO 

A)i 


(1 + R) 2 e~ 4R 
R 


5[« 


-2B(R-r) _ 2e -B(R-r) 




-2B(R-r) -B(R- r 


- e 


>dR 


(17a) 
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respectively, where z = 1, a Q = 1 (using atomic units for the case of H(ls)), 
where from equations (3), (14), and (15) 


-o 1 n a 
R 03 " p ln e 


and 


form 


-I ln ( 1+ V*l) 

The integration indicated in equations (16a) and (17a) may he put in the 


V = T 


g 3 = h(4,R 03 ) - | h(4 + p,Ro 3 ) + k(4+ P,^ 3 ) 


(16ns) 


and 


Si = h(4,%) - | e 2Br h(4 + 2B,R 01 ) + 2 | e Br h(4 + B,^) 

+ | Be 2Br k(4 + 2B,R 01 ) - 2 Be Br k(4 + B,R Q1 ) (17b) 


where h and k are defined hy 


h(s,t) = 


J /*oo 

: ^ 


2 -sX 
X) e 


dX = -Ei(-st) + - e“ st + ~ e“ st (st + l) 
s 


and 


/ 


k(s,t) = / (1 + X) 2 e" sX dX = e' 


■st 


a. * i)z + 2 ii+ii + 4 

s s 2 s 3 


where -Ei(-X) is the exponential integral tabulated hy Jahnke and Emde 
(ref. 7). 

Ionization Cross Section 


The ionization cross section for the interaction along a triplet po- 

tential and Oj. for the interaction along a singlet potential are given by- 
equations (12), (16b), and (17b). These cross sections are presented in fig- 
ure 3. Also shown in figure 3 is the cross section "5 which would be observed 
experimentally if there were no mechanism for systematic selection of spin 
states. This average cross section *d is given by (ref. 8) 


7 



4 


( 18 ) 



the incident particle at the large values 
the approximation R c ^ b c can be used. 


1 3 

a = — a, + -T a* 
4 1 4 3 


CLASSICAL PATH APPROXIMATION 

The classical collision of 
particles under the influence of 
certain central potentials (those 
which decrease faster than l/r as 
r -> oo ) follow paths for which the 
angular deviation decreases as im- 
pact parameter increases. At some 
critical angle classical mechanics 
is no longer valid (ref. 9), and 
there exists a maximum impact param- 
eter b c , beyond which quantum ef- 
fects on "particle path" must be 
considered. Rosen 1 s approximation 
requires that the particle paths be 
treated classically; hence, for 
this approximation to be valid all 
inelastic transitions must occur at 
distances of closest approach less 
than R , where R c corresponds to 
the impact parameter b c . Since 
there is very little deflection of 
of impact parameter where b c occurs, 


The relation between R c and the collision parameters for an exponential 
potential of the form 


-Cpr 

V(r) = C-^e 


( 19 ) 


may be obtained from Mason, Vanderslice, and Raw (ref. 10 ). Their equations 
(l) and (69) combine to give 


- | ln ( C 2 B c> + C 2 E c 



( 20 ) 


where k = pv/fi is the wave number of the incident atom. 

•Z 

For the case of H(ls) - H(ls) collisions following the Z u interaction 
potential (eq. (14)) R c is about 3.5 a Q . The 1 Z g interaction potential at 
long range can be approximated by an exponential attractive potential (eq. (19) 
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Cg = B). This procedure gives a value of about 5a Q 


-with CL = -2De B ^, 

for R7 
c 

The effective 


range of the interaction may he approximated by 



( 21 ) 


where R eff < la for e < 21 (fig. 3). Hence R c > R eff for the H(ls) 
collisions reported herein. 


CONCLUDING REMARKS 

It is surprising that more work has not been done with Rosen's approxima- 
tion, since it applies to the interesting thermal plasma energy regime just 
above the ionization threshold. This is an energy regime which is difficult to 
investigate experimentally. 

The merit of the method may be judged from a comparison made in references 
1 and 2 with the experimental data of Rostagni (ref. 11) for the case of helium- 
helium and argon-argon collisions. While the experiments were performed at col- 
lision energies where the validity of the sudden approximation is questionable, 
the theoretical and experimental cross sections agreed to within a factor of 
about two. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, September 22, 1965. 
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APPENDIX - SYMBOLS 


An 


distance from electron 1 to nucleus A 


A 2 

a 

a o 

B 

B 1 

B 2 

b 

°i,c 2 

D 

e 

f 

g l’ g 3 

h(s,t) 

ft 

I 


I l> I 2 > I 3’ 

h’h 


distance from electron 2 to nucleus A 


k(s, t) 


m 


atomic parameter, — 


ft 


r 2ml 


radius of first Bohr orbit, 0.529 A 

constant in Morse potential 

distance from electron 1 to nucleus B 

distance from electron 2 to nucleus B 

impact parameter 

critical impact parameter 

constants in exponential potential 

dissociation energy in Morse potential 

electron charge 

integral defined by eq. (10) 

integrals defined by eqs. (16a) and (17a) 


/ 


(1 + X)V sX 

X 


ax 


Planck's constant divided by 2jt 
ionization energy 

Rosen's integrals evaluating eq. ( 8 ) 


wave number, 

7 ft 


i: 


(1 + x) 2 e" sX ax 


electron mass 

ejected electron momentum 
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R 


R eff 

R or r o3 


r 


r 


12 


s 


t 

Uj/R) 

V R > 

V 


V(r) 

V R > 


V 


X 




a 

P 

r 

e 

M- 

a 


interatomic distance of closest approach 
critical distance of closest approach 
effective range of inelastic transition 

distance of closest approach corresponding to b = 0 for interaction 
potentials U^(R) and U 3 (R), respectively 

radial coordinate 

distance from electron 1 to electron 2 
variable in integral functions h(s, t), k(s, t) 
variable in integral functions h(s, t), k(s, t) 
singlet hydrogen interaction potential 
triplet hydrogen interaction potential 
perturbation potential (eq. (5)) 
exponential potential 

perturbation matrix element defined by eq. (4) 
velocity of incident atom 
dummy variable of integration 
ionization parameter, y/i 1 
effective charge 

range constant in triplet potential 
energy constant in triplet potential 

radius of energy minimum in Morse potential (singlet potential) 

1 9 

collision energy, — \iv 
reduced mass of colliding atoms 
ionization cross section 


constant, 


64 g 4 (t) k£ 

(2ml) 5 / 2 ! 2 
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singlet ionization cross section 
triplet ionization cross section 
a average ionization cross section defined by eq. ( 18 ) 

t variables of integration representing volume 

cp(r) normalized hydrogen Is wave function 

\|/ f wave function for final state of colliding atoms given by product of Is 

wave function and plane wave representing ejected electron 

^ wave function for initial state of colliding atoms given by product of 
two Is wave functions 


12 



REFERENCES 


1. Rosen, Philip: Low-Energy Inelastic Atomic Collisions. Phys. Rev., vol. 

109, no. 2, Jan. 15, 1958, pp. 348-350. 

2. Rosen, Philip: Ionization Cross Section of Argon-Argon Collisions Rear 

Threshold. Phys. Rev., vol. 109, no. 2, Jan. 15, 1958, pp. 351-355. 

3. Bates, D. R. ; and Griffing, G. : Inelastic Collisions Between Heayy 

Particles. I. Excitation and Ionization of Hydrogen Atoms in Fast 
Encounters with Protons and with Other Hydrogen Atoms. Proc. Phys. Soc., 
sec. B, vol. 66, Nov. 1953, pp. 961-971. 

4. Petit Bois, G. : Tables of Indefinite Integrals. Dover Pub., Inc., 1961. 

5. Slater, John C.: Quantum Theory of Molecules and Solids. Vol. I. Elec- 

tronic Structure of Molecules. McGraw-Hill Book Co., Inc., 1963. 

6. Dalgarno, A.; and Lynn, N. : Resonance Forces at Large Separations. Proc. 

Phys. Soc., sec. A, vol. 69, Nov. 1956, pp. 821-829. 

7. Jahnke, E.; and Emde, F. : Tables of Functions. Fourth ed., Dover Pub., 

Inc., 1945. 

8. Konowalow, Daniel D.; Hirschfelder, J. 0.; and Linder, Bruno: Low- 

Temperature, Low-Pressure Transfer Coefficients for Gaseous Oxygen and 
Sulfur Atoms. J. Chem. Phys., vol. 31, no. 6, Dec. 1959, pp. 1575-1579. 

9. Burhop, E. H. S. : Theory of Collisions. Vol. 1 of Quantum Theory, D. R. 

Bates, ed. , Academic Press, Inc., 1961, ch. 9. 

10. Mason, E. A.; Vanderslice, J. T. ; and Raw, C. J. G. : Quantum Effects in 

Small-Angle Molecular-Beam Scattering. J. Chem. Phys., vol. 40, no. 8, 
Apr. 15, 1964, pp. 2153-2164. 

11. Rostagni, A.: Positive and Neutral Rays. III. Ionisation by Collision 

of Atoms. Ruovo Cimento, vol. 11, Nov. 1934, pp. 621-634. 


NASA -Langley, 1965 


E-3106 


13 



"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 

important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington D.C. 20546 



